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Spatially  resolved  rotational  temperatures  have  been  obtained  within  the  boundary  layer 
of  a  flat  plate  model  in  a  Mach  5.1  flow  using  emission  spectroscopy.  The  temperatures  were 
obtained  by  matching  the  measured  nitrogen  second  positive  0-2  rovibrational  band  with  a 
calculated  one.  Temperature  profiles  are  given  above  the  cathode  and  the  anode.  The 
maximum  temperature  obtained  above  the  cathode  did  not  correspond  to  the  surface  of  the 
model,  but  rather  at  an  elevation  0.55  mm  above  the  surface.  This  indicates  that  heat  is 
traveling  from  the  discharge  into  the  plate  for  a  period  of  time  after  the  discharge  is  ignited. 
This  characteristic  in  the  temperature  profile  sheds  light  on  some  of  the  effects  a  plasma  has 
on  the  flow  field  as  determined  by  pitot  probe  measurements  and  total  lift  measurements. 
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I.  Introduction 

POSSIBLE  benefits  of  magneto-aerodynamics,  which  are  of  great  interest  to  the  aerospace  community  include 
plasma  flow  control  with  applications  to  boundary  layer  transition  control,1  shock  boundary  layer  interaction 

control,  "  3  virtual  flap  actuators  for  aerodynamic  control,  4  drag  reduction  in  flight, 5  6  and  steering  of  high-speed 

7 

aircraft.  These  benefits  may  be  achievable  by  generating  a  plasma  discharge  and  applying  a  magnetic  field  in  high¬ 
speed  flows;  however,  they  cannot  be  realized  fully  until  the  properties  of  plasma  located  in  high-speed  flows  are 
well  characterized.  At  the  present  time  the  combination  of  both  Joule  heating  effects  and  heating  from  the  electrode 
surface  through  convection  have  shown  to  noticeably  alter  the  bulk  of  the  fluid  flow  and  alterations  due  to  a 

Lorentzian  force  are  not  discemable.  This  has  led  to  the  work  being  presented  in  this  paper  which  is  the 
characterization  of  the  temperature  profiles  present  in  the  boundary  layer.  This  is  done  by  spatially  resolving 
rotational  temperatures  at  different  positions  within  the  boundary  layer  using  emission  spectroscopy. 

Hypersonic  flow  characteristics  above  the  flat  plate  change  significantly  due  to  heating  caused  from  igniting  a 
plasma.  These  changes  are  understood  by  measuring  the  thermal  boundary  layer.  There  are  a  limited  number  of 
techniques  capable  of  performing  these  measurements  in  nonequilibrium  plasmas.  One  technique  is  to  insert  a 
temperature  probe  into  the  discharge.  This  method,  which  measures  the  recovery  temperature,  is  intrusive  and  will 
not  provide  accurate  measurements  if  it  provides  results  at  all  because  a  probe  such  as  a  thermocouple  interacts  with 
the  self  induced  electric  field  of  the  discharge.  Also  temperatures  measured  with  a  thermocouple  will  be  skewed 
hotter  than  what  is  physically  possible  due  to  electrons  condensing  on  the  probe  surface.  These  limitations  with  an 
intrusive  temperature  measurement  demonstrate  the  need  for  a  nonintmsive  technique,  such  as  spectroscopy.  A 
spectroscopic  method  used  for  many  nonequilibrium  plasmas  is  rotational  spectroscopy.  Many  times  the  rotational 
temperature  and  the  translational  temperature  in  a  plasma  are  close  to  being  equal.  However,  to  measure  rotational 
temperatures  spectroscopically  the  rotational  modes  do  not  have  to  be  in  equilibrium  with  the  translational, 
vibrational  and  electronic  modes,  the  rotational  modes  only  need  to  be  in  self  equilibrium.  If  this  is  true  then  the 
rotational  temperature  can  be  found  for  complicated  band  structures  by  matching  a  calculated  spectrum  with  the 
measured  spectrum. 10-12  Another  technique  is  the  Boltzmann  plot  method. I3' 14 

II.  Test  Facilities  and  Equipment 

Tests  were  conducted  within  a  blow-down,  free  jet,  Mach  5.1  flow  channel  located  at  Wright  Patterson  Air  Force 
Base  (see  Fig.  1).  The  flow  channel  is  constructed  from  acrylic  plastic  and  assembled  with  nylon  screws.  The  test 
section  has  a  removable  quartz  cover  plate  to  allow  transmittance  of  radiation  in  the  near  infrared  and  visible 
wavelengths.  The  stagnation  temperature  is  weakly  dependent  on  the  outdoor  air  temperature  and  is  around  270K. 
The  stagnation  pressure  has  an  operating  range  of  0. 1  to  1.0  atm.  More  detail  on  the  flow  channel  characteristics  can 
be  found  in  a  prior  publication. 15 

A  flat  plate  model  made  out  of  a  dielectric  material  with  copper  electrodes  was  used  in  this  study.  The  cathode  of 
the  model  was  positioned  up  stream  for  every  run.  A  schematic  of  the  model  is  shown  in  Fig.  2. 

The  emitted  radiation  from  the  surface  of  the  model  is  collected  by  a  Newport  KBX157  bi-convex  lens  with  an 
effective  focal  length  of  125  mm.  The  lens  is  positioned  two  times  the  focal  length  away  from  the  discharge.  A 
Newport  aperture  stop  is  used  in  front  of  the  lens  to  help  eliminate  stray  radiation.  The  image  of  the  discharge  is 
focused  at  a  distance  of  two  times  the  focal  length  of  the  collection  lens  onto  a  Roper  Scientific  LG-455-020-3  fiber 
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Figure  2.  Model  schematic.  All  dimensions  are  in  cm. 


2 

optic  bundle.  The  fiber  optic  bundle  channels  collected  radiation  from  a  spot  size  of  0.125  mm  into  an  Acton  2756 
spectrometer  where  it  is  diffracted  and  focused  onto  an  Andor  DU  440-BU  CCD  camera.  The  spatial  resolution  of 
the  measuring  technique  is  given  by  the  spot  size  divided  by  the  magnification  which  is  one  for  the  lens  used.  The 
spatial  resolution  is  0. 125  mm.2  The  spectrometer  has  a  spectral  resolution  of  0.013  nm  at  380.415  nm  using  a  3600 
grooves  per  mm  holographic  grating  blazed  at  250  nm.  The  slit  height  and  width  respectively  are  4  mm  and  13pm. 
The  optical  components  are  mounted  onto  a  rail  positioned  perpendicular  (x-direction)  to  the  flow  through  the  wind 
tunnel.  This  rail  is  mounted  onto  another  rail  aligned  in  the  flow  direction  (y-direction).  This  allows  measurements 
to  be  taken  at  different  positions  in  the  y-direction.  A  scale  is  scribed  onto  the  y-direction  rail  so  that  accurate 
positioning  can  be  done  in  the  flow  direction.  The  y-direction  rail  is  mounted  to  a  steel  elevating  table. 
Measurements  at  different  elevations  in  the  discharge  are  obtained  by  raising  and  lowering  the  table.  The  table 
elevates  by  turning  a  wheel  with  diameter  of  61.3  cm.  Each  complete  turn  of  the  wheel  corresponds  to  a  change  in 
elevation  of  8.3  mm.  The  wheel  is  divided  into  61  even  divisions  where  each  division  corresponds  to  a  change  in 
elevation  of  0.14  mm.  Fig.  3  shows  the  table  along  with  the  optics  and  how  they  are  oriented  with  the  test  section. 
Alignment  of  the  table  with  respect  to  the  model  was  done  using  a  large  square.  A  level  was  used  to  adjust  the  table 
and  model  so  that  these  surfaces  are  parallel. 
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Figure  3.  Photograph  of  the  optical  system,  the  table  and  the  orientation  of  each  piece  with  respect  to  the  test 

section. 


ITT.  Rotational  Temperature  Determination 

The  rotational  temperature  was  obtained  by  comparing  the  measured  spectrum  of  the  nitrogen  second  positive  0- 
2  transition  to  different  theoretical  spectrums  calculated  for  different  rotational  temperatures.10"1’  The  rotational 
temperature  obtained  corresponds  to  the  calculated  spectrum  that  best  matches  the  experimental  spectrum.  This 
procedure  is  performed  by  a  program  called  N2SPECF1T. 1 1 

Assuming  an  optically  thin  collection  cone  is  used;  the  intensity  of  a  rovibrational  line  can  be  written  as16 
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The  approximation  to  the  measured  spectrum  is  obtained  from  Eqn.  (1)  by  summing  up  all  the  rovibrational  lines 
and  convolving  them  with  the  instrumental  line  shape  function  P(X,  A') .  The  different  broadening  mechanisms  of 
the  measured  spectrum  produce  half-widths  that  are  at  least  an  order  of  magnitude  smaller  than  the  half-width  of  the 
instrumental  line  function  in  the  environment  being  investigated.  Therefore  the  half-widths  of  the  different 
broadening  mechanisms  present  in  the  measured  spectrum  can  be  approximated  by  a  Dirac  delta  function  with 
minimal  error.  This  assumption  makes  the  convolution  easier.  The  instrumental  line  shape  function  was  obtained  by 
measuring  and  curve  fitting  the  spectra  of  a  mercury  discharge  tube  at  435.833  nm.  This  value,  when  convolved, 
may  not  broaden  the  calculated  spectra  enough  to  match  the  measured  due  to  the  assumption  made  with  the 
broadening  mechanisms.  The  variance  of  the  calculated  spectrum  with  respect  to  the  measured  spectrum  can  be 
improved  by  modifying  the  instrumental  line  function  obtained.  This  produces  negligible  error.  The  calculated 
spectrum  is  then  given  by 
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The  method  described  above  to  measure  rotational  temperatures  assumes  that  the  rotational  temperature  is 
uniform  in  the  collection  volume.  This  is  clearly  not  the  case  in  the  current  paper  since  the  measurements  were  taken 
in  a  spatially  non-uniform  plasma.  This  error  was  reduced  by  minimizing  the  size  of  the  collection  volume.  Also 
numerical  verification  was  used  to  understand  what  temperature  value  would  be  obtained  from  a  collection  volume 
described  by  a  range  of  rotational  temperatures.  This  was  done  by  generating  different  numerical  spectra  that 
contained  equal  weighted  contributions  from  two  temperatures  and  then  using  these  spectra  as  inputs  to  the  program 
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N2SPECFIT.  The  results  are  shown  in  Table  1  for  seven  different  combinations  of  two  temperature  spectrums 
separated  by  50  K.  The  tabulated  values  indicate  that  the  software  converges  on  a  value  around  the  average 
temperature  of  the  combined  temperature  spectra.  Therefore  the  rotational  temperatures  measured  in  this  paper 
should  be  close  to  an  average  temperature  over  the  spot  size  the  radiation  was  collected  from. 


Temperature  Mixture  (K) 

Converged 

Temperature  Value  (K) 

450  and  400 

423.9 

400  and  350 

373.9 

350  and  300 

323.8 

300  and  250 

273.7 

250  and  200 

223.5 

200  and  150 

172.8 

150  and  100 

120.7 

100  and  50 

78.5 

Table  1.  Converged  temperatures  found  for  different  spectra  mixtures. 


IV.  Experimental  Procedures 


A.  Calibration 

The  wavelength  falling  onto  the  center  pixel  of  the  Andor  CCD  camera  was  calibrated  using  a  helium  discharge 
tube.  The  calibration  of  the  Acton  spectrometer  assumes  the  wavelength  focused  onto  the  center  pixel  is  dispersed 
linearly  with  the  rotating  grating.  The  software  that  operates  the  spectrometer  takes  the  peak  position  of  two 
different,  single  electronic  transitions  to  calculate  the  offset  and  slope  used  for  the  calibration.  The  y-offset  for  the 
calibration  was  found  using  the  helium  line  at  318.774  nm  and  the  slope  using  the  helium  line  at  447.147  nm. 
Several  other  electronic  transitions  were  used  to  verify  the  calibration.  The  slope  represents  the  rate  of  change  of  the 
wavelength  over  the  angle  of  the  grating.  The  calibration  was  off  by  five  pixels  in  the  spectral  region  of  interest  so 
the  offset  value  was  adjusted  by  this  amount.  The  wavelength  dispersion  across  the  other  pixels  of  the  CCD  camera 
was  found  using  the  grating  equations.  These  equations  assign  a  wavelength  to  each  pixel  and  are  dependent  on  the 
geometry  of  the  spectrometer.  The  grating  equations  are  given  by 
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The  response  of  the  camera  to  a  constant  intensity  source  was  not  uniform  across  the  CCD  and  was  calibrated 
using  a  tungsten  filament  lamp  with  known  spectral  emission.  This  is  a  crucial  step  for  properly  measuring  the 
intensity  ratios  of  the  different  rotational  transitions  of  a  band  which  is  necessary  for  obtaining  accurate  temperature 
measurements.  Signal  due  to  stray  light  inside  the  spectrometer  caused  from  the  imperfections  of  the  mirrors  used 
was  accounted  for  by  using  a  filter  that  blocks  the  wavelength  range  over  which  the  calibration  curve  is  being 
obtained.  The  calibration  curve  is  given  by 


Calibration  Curve(X)  = 


Tungsten  Lamp  Spectra(X)~  Stray  Light  Spectra(X) 
Theoretical  Tungsten  Lamp  Spectra(X) 


(6) 
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The  calibrated  data  is  given  by 


Calibrated  Data{X)  = 


Raw  Data{/ l) 
Calibration  Curve(A ) 


(7) 


B.  Alignment  and  Preparation 

Careful  alignment  of  the  table  with  respect  to  the  model  was  performed.  A  large  square  was  used  to  orient  the 
table  perpendicular  to  the  model.  A  level  was  used  to  make  sure  both  the  model  and  table  surfaces  were  horizontal. 
This  is  crucial  and  it  was  found  by  calculation  that  half  a  degree  difference  leads  to  significant  error.  The 
measurement  volume  of  the  optics  inside  the  test  section  was  found  by  shinning  a  laser  through  the  back  side  of  the 
fiber  optical  bundle.  This  technique  was  also  used  to  locate  the  position  of  the  leading  edge  of  the  model  relative  to 
the  y-direction  rail  and  this  was  taken  as  the  position  where  y=0.  The  fiber  optic  cable  was  then  carefully  set  back 
into  its  original  position  in  the  entrance  slit  of  the  spectrometer. 

The  camera  was  cooled  to  218  K  to  reduce  dark  signal  noise.  All  spectra  scans  were  acquired  with  full  vertical 
binning,  16  second  pixel  read  out  time,  and  background  corrected  spectra.  The  counts  for  each  scan  were  held  at  a 
constant  value  just  under  saturation  by  increasing  the  exposure  time.  This  means  the  overall  exposure  time  was 
different  for  points  measured  at  different  positions  in  the  discharge.  Each  scan  was  accumulated  four  times  to 
increase  the  signal  to  noise  ratio. 


V.  Results 

Rotational  temperatures  were  taken  over  the  anode  at  a  current  of  25mA  and  over  the  cathode  at  25mA  and 
50mA.  The  stagnation  pressure  was  held  at  370  torr  for  all  cases  and  the  stagnation  temperatures  varied  between 
274  K  and  281  K.  The  profdes  over  the  cathode  are  shown  in  Fig.  4.  Points  below  an  elevation  of  0.3  mm  include 
part  of  the  model  surface  in  the  optical  line  of  sight  and  the  temperatures  are  skewed  colder,  possibly  due  to  some 
reflected  radiation  off  of  the  model  surface.  These  points  were  omitted  from  the  profdes  and  the  first  point  taken 
corresponds  to  an  optical  line  of  sight  free  of  the  model  surface.  The  first  point  of  the  50mA  case  has  a  temperature 
of  360  K  at  an  elevation  of  0.3  mm  and  increases  to  a  maximum  value  of  384  K  at  an  elevation  of  0.55  mm.  The 
negative  slope  in  the  temperature  profile  at  the  model  surface  shows  that  heat  is  traveling  into  the  cathode.  The 

dT 

slopes  -  at  the  cathode  surface  for  the  25mA  and  50mA  cases  are  108  and  165  K/mm  respectively.  Assuming  the 

dy 

rotational  temperature  is  in  equilibrium  with  the  translational  temperature  this  is  an  important  observation  because  it 
is  a  possible  explanation  for  some  transient  behavior  seen  in  pitot  probe  measurements17  and  total  lift 
measurements18  made  in  this  facility.  The  assumption  that  translational  energy  modes  are  in  equilibrium  with  the 
rotational  modes  stems  from  the  small  energy  separation  of  a  rotational  transition  which  allows  for  rapid  energy 
exchange  between  the  two  energy  storage  modes.19 

Outside  of  the  thermal  boundary  layer  the  temperature  is  measured  to  be  82  K.  This  temperature  is  higher  than 
the  analytically  determined  free  stream  static  temperature  downstream  of  the  shock.  The  source  of  this  error  is  not 
fully  understood,  but  is  believed  to  be  caused  by  a  combination  of  a  weak  signal  in  this  region  of  the  discharge  and 
reflected  radiation.  The  reflected  radiation  is  believed  to  be  coming  from  hotter  portions  of  the  discharge,  thus 
raising  the  measured  temperature. 

The  final  interesting  feature  of  the  temperature  profile  is  the  change  recorded  across  the  shock.  This  change  can 
be  seen  in  Fig.  4  but  it  is  easier  to  see  in  Fig.  5  which  is  an  enlarged  view  of  the  shock  region.  The  temperature  drop 

calculated  for  an  oblique  shock  with  shock  angle  of  18.3  °  and  a  Mach  number  of  5.1  is  18  K.  The  shock  angle  was 
obtained  by  using  a  protractor  to  measure  the  angle  from  a  Schlieren  image.  The  measured  value  found  is  only  7  K. 
The  discrepancy  is  not  fully  understood  but  falls  under  the  same  explanation  given  for  the  skewed  temperature 
values  recorded  in  the  free  stream.  The  temperature  jump  across  the  shock  is  shown  as  having  occurred  over  a 

2 

1.5  mm  distance.  This  discrepancy  is  from  both  the  spatial  resolution  of  0.125  mm  and  the  accuracy  of  how  level 
the  optical  system  is  to  the  model.  Fig.  6  shows  one  of  the  measured  radiation  spectmms  from  this  region  overlaid 
with  the  calculated  spectrum  from  N2SPECFIT.  The  calculated  spectrum  appears  to  agree  well  with  the  measured 
spectrum. 
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The  rotational  temperature  profile  for  the  anode  is  shown  in  Fig.  7.  The  glow  over  the  anode  is  a  lot  smaller  than 
that  for  the  cathode.  Because  the  discharge  over  the  anode  only  emits  significant  radiation  for  a  small  height  above 
the  anode  surface,  measurements  of  the  rotational  temperature  could  only  be  made  up  to  a  location  3  mm  above  the 
surface.  The  portion  of  the  temperature  profile  which  could  be  obtained  shows  a  constant  temperature  of  about 
299  K. 

It  is  ideal  to  measure  a  single  plasma  temperature  instead  of  an  averaged  value  over  a  small  collection  volume. 
However,  in  our  system  there  are  three  means  by  which  a  temperature  gradient  is  present  in  the  measurement.  As 
mentioned  before  there  is  a  gradient  present  at  different  elevations  above  the  plate  surface.  The  results  in  Fig.  4 
show  this.  This  occurs  because  a  small  but  finite  measurement  size  in  the  vertical  direction  is  used.  Another  gradient 
exists  in  the  spanwise  direction  of  the  discharge.  The  third  temperature  gradient  is  detected  when  the  air  heats  up 
from  the  discharge  over  time  during  the  accumulation  process.  This  occurs  from  the  electrodes  heating  up  over  time. 
The  accumulation  time  used  in  this  work  varied  depending  on  the  location  in  the  discharge.  For  locations  close  to 
the  model  surface  short  accumulation  times  varying  from  6  to  12  seconds  were  used  because  the  light  intensity  was 
strong.  For  locations  away  from  the  plate  longer  accumulation  times  up  to  90  seconds  were  used  because  the 
intensity  for  the  light  signal  was  weak.  From  the  results  shown  in  Figs.  8  and  9  it  can  be  seen  that  the  large  changes 
in  temperatures  with  time  occur  close  to  the  plate  surface.  Away  from  the  surface  the  changes  are  much  smaller. 
Figure  8  shows  the  results  for  temperature  measurements  taken  over  time  in  the  free  stream  at  an  elevation  of  4  mm. 
The  maximum  change  in  temperature  is  2  K  and  corresponds  to  50  mA  of  current  with  the  discharge  on  for  40 
seconds.  The  last  point  taken  in  Fig.  9  for  the  50  mA  case  shows  the  temperature  decreasing  instead  of  increasing. 
This  is  due  to  the  discharge  breaking  down.  Figure  9  shows  the  temperature  increasing  as  a  function  of  time  in  the 
boundary  layer.  Notice  that  the  temperature  increase  in  the  boundary  layer  is  small  for  short  time  durations. 
Consequently  the  exposure  times  in  the  boundary  layer  where  the  signal  is  very  strong  were  small  and  in  the  free 
stream  where  the  signal  was  weak  were  long.  Thus  it  is  concluded  that  the  variation  of  temperature  with  time  is  an 
unimportant  effect  in  these  measurements. 


Figure  4.  Rotational  temperature  profiles  over  the  cathode  at  currents  of  25  mA  and  50  mA. 
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Figure  5.  Close  up  of  temperature  jump  across  the  shock  from  Fig.  4. 
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Figure  6.  Output  from  N2SPECFIT  shows  a  measure  rovibrational  band  overlaid  with  the  calculated  fit. 
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Figure  8.  Time  dependence  on  temperature  in  the  free  stream  at  the  cathode  for  an  elevation  of  4.09  mm. 
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Figure  9.  Time  dependence  on  temperature  in  the  boundary  layer  at  the  cathode  for  an 

elevation  of  0.273  mm. 


VI.  Conclusion 

Spatially  resolved  rotational  temperatures  have  been  obtained  inside  the  boundary  layer  of  a  flat  plate  model 
with  a  surface  plasma  discharge  in  a  Mach  5.1  air  flow.  The  maximum  temperature  measured  for  both  the  25mA 
case  and  50mA  case  occurs  over  the  cathode  at  an  elevation  of  0.546  mm.  The  maximum  temperatures  for  these  two 
different  currents  were  384  K  and  301  K  respectively.  Because  the  temperature  drops  in  both  vertical  directions 
from  this  maximum  it  can  be  conclude  that  heat  is  being  convected  into  the  surface  of  the  plate  for  a  period  of  time 
when  the  discharge  is  first  turned  on.  This  is  a  possible  explanation  for  some  trends  seen  in  pitot  and  total  lift 
measurements  made  in  other  work  done  by  the  authors.  The  temperature  profile  over  the  cathode  shows  a  small 
temperature  change  across  an  oblique  shock.  The  measured  temperature  change  across  the  shock  of  7  K  does  not 
agree  with  the  calculated  change  of  18  K.  The  disagreement  is  believed  to  be  caused  by  a  weak  signal  in  this  region 
of  the  discharge  and  from  reflected  light.  The  temperature  profile  obtained  for  the  anode  shows  a  relatively  constant 
temperature  value  of  299  K  out  to  an  elevation  of  3  mm.  This  work  also  shows  that  the  temperature  of  the  gas  above 
the  cathode  does  change  with  time.  This  change  with  time  diminishes  at  locations  well  above  the  copper  electrode. 
This  temporal  dependence  along  with  the  negative  temperature  slope  at  the  cathode  surface  indicates  that  the  copper 
cathode  is  heating  and  thus  reducing  the  amount  of  energy  being  dumped  into  the  air  flow,  for  a  period  of  time,  right 
above  the  model  surface. 
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